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FOAM  GENERATION  AND  AIR  ENTRAINMENT 
NEAR  A  FREE  SURFACE 


INTRODUCTION 

The  flow  in  a  breaking  deep  water  wave  is  a  type  of  agitation  which  is  conducive  to  foam  for¬ 
mation.  The  downstream  wake  of  a  surface  vessel  moving  in  the  seaway  also  is  characterized  by  a 
layer  of  foamy,  agitated  water  near  the  surface,  and  a  subsurface  layer  of  bubbles  or  a  bubble  cloud. 
An  extensive  study  of  ship  wakes  and  white  water  dating  from  the  years  during  and  just  after  World 
War  Two  is  given  in  a  U.S.  Navy  report  compiled  by  the  National  Defense  Research  Council  (1969). 
More  recently  two  studies  have  been  conducted  at  NRL  to  define  more  clearly  the  white-water  region 
of  a  surface  ship  wake  (Peltzer,  1984;  Peltzer,  Garrett  and  Smith,  1987).  The  highly  agitated  region 
of  white  water  just  downstream  from  a  typical  surface  ship  is  shown  in  Fig.  1. 

The  purpose  of  this  report  is  to  discuss  the  basic  features  of  foam  and  bubble  formation,  and  air 
entrainment  as  these  features  are  known  at  the  present  time.  We  shall  discuss  possible  means  for  relat¬ 
ing  these  basic  features  to  the  processes  of  deep  water  wave  breaking  and  of  wake  formation  and 
development  for  surface  ships  moving  in  the  seaway. 

PROPERTIES  OF  FOAM  AND  BUBBLES 

A  foam  is  a  coarse  dispersion  of  gas  in  a  liquid,  with  most  of  the  volume  being  gas.  The  liquid 
is  in  thin  sheets  or  lamellae  between  gas  bubbles.  Foams  are  formed  by  dispersing  gas  into  a  liquid 
either  by  reducing  the  pressure  over  the  liquid  containing  a  dissolved  gas,  or  by  mechanical  gas-liquid 
contact.  Mechanical  contact  can  be  accomplished  by  pumping  gas  through  a  porous  septum  or  a  capil¬ 
lary,  or  by  strongly  agitating  the  liquid.  There  must  be  some  form  of  stabilizer  or  surfactant  material 
present  if  the  foam  is  not  to  break  or  collapse  immediately.  This  means  that  foams  cannot  be  formed 
from  pure  liquids.  The  foamy  region  of  a  highly  agitated  breaking  wave  or  surface  wake  region  is 
confined  to  the  near-surface  bubbly  region.  The  remainder  of  the  bubble  distribution  or  cloud  extends 
to  greater  depths  due  to  the  turbulent  motion  of  the  water. 

DeVries  (1972)  describes  the  structure  of  foams.  If  the  foam  contains  spherical  bubbles,  up  to 
74  percent  of  the  foam  can  be  gas  if  the  bubbles  have  a  distribution  of  sizes.  The  surface  area  A  of 
the  gas-liquid  interface  is  given  by 
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Fig.  1  —  A  typical  aerial  view  of  the  white  water  region 
behind  a  surface  ship  underway;  from  Peltzer  (1984) 
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The  factor  (£  n ;d{3)  is  the  reciprocal  of  the  mean  bubble  diameter,  which  depends  on  the 

bubble  size  distribution.  As  a  foam  ages,  the  liquid  drains  from  between  the  bubbles,  and  the  bubbles 
begin  to  coalesce.  This  causes  the  surface  area,  A,  to  decrease  with  time. 

As  the  liquid  drains,  the  bubble  walls  become  thin,  and  the  bubbles  become  polyhedral.  Then 
the  lamellae  separating  them  are  distorted  because  of  surface  tension  effects.  Pressure  differences  exist 
across  strongly  curved  interfaces,  so  that 


Ap 


R{  R2 


(2) 


where  7  =  surface  tension 

Ri,  R2  =principal  radii  of  curvature  of  the  interface. 

There  are  illustrations  of  the  meaning  of  the  principal  radii  in  Fig.  2  for  two  isolated  spherical  bub¬ 
bles,  and  in  Fig.  3  for  two  bubbles  which  have  a  very  thin  lamella  between  them.  It  is  observed  that 
in  foams  having  polyhedral  bubbles,  the  bubble  walls  meet  in  threes,  and  each  angle  is  120.  This 
angle  does  not  depend  on  bubble  size  or  on  bubble  size  distribution.  The  reason  suggested  is  that  the 
pressure  difference  between  bubbles  is  determined  mainly  by  the  smallest  radius  of  curvature. 


The  liquid  in  a  foam  is  located  in  the  so-called  plateau  borders  where  the  walls  of  the  polyhedral 
bubbles  join.  This  is  illustrated  in  Fig.  4.  Gravity  forces  and  capillary  suction  can  produce  very  com¬ 
plicated  flow  patterns  in  the  bubble  walls  as  liquid  drains  from  the  foam.  The  cross-sectional  area  for 
flow  through  the  plateau  borders  in  a  foam  having  dodecahedral  bubbles  is 
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The  notation  here  is  that  used  in  Fig.  4.  The  gas-liquid  ratio  for  this  situation  is 

Vr  60A pb  E  n.d;  E  n:d? 

~ —  =  - -  - T  30  - 7.  (4) 

VG  2.445  t  £  md>  E  n,43 

Here  dt  is  the  diameter  of  a  spherical  bubble  before  it  becomes  a  dodecahedron. 

The  bubble  size  distribution  is  an  important  factor  used  in  describing  the  structure  of  foams. 
This  can  be  measured  experimentally  by  determining  the  frequency  distribution  of  bubble  sizes  photo¬ 
graphically  in  a  plane  passing  through  the  foam.  To  produce  such  a  plane,  the  foam  can  be  placed  in 
a  container  having  a  flat,  transparent  wall  or  the  foam  can  be  frozen  and  part  of  it  sliced  away.  It  has 
been  shown  that  the  size  distribution  is  the  same  in  any  plane  in  a  foam,  because  the  bubbles  are  ran¬ 
domly  distributed.  The  frequency  distribution  of  bubble  sizes  observed  in  a  plane  will  differ  by  a  con¬ 
stant  from  the  frequency  distribution  over  a  certain  foam  volume.  This  is  discussed  by  deVries  and 
Lemlich  (1972).  The  two  frequency  distribution  functions  are  related  by 


nsf(r)dr  =  2  rCsNF(r)dr, 


(5) 


where 
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=  number  of  bubbles  per  unit  area  of  the  plane; 

=  number  of  bubbles  per  unit  volume  in  the  foam; 
=  size  distribution  function  for  the  plane; 

=  size  distribution  function  for  the  foam  volume; 

=  bubble  radius. 
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Fig.  3  —  Interference  geometry  for  bubbles  close  together  with  a  thin  lamella 
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Integrating  between  r  =  0  and  r  =  oo  relates  n  s  and  N  to  the  mean  bubble  radius  <  r  > , 

ns=2CsN<r>.  (6) 


The  frequency  distributions  are  related  by 


fir)  —  F(r). 

<r> 


(7) 


The  constant  Cs  can  be  related  to  the  mean  bubble  volume,  (4t/3)  <  r3  > ,  by 
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or 


C.  = 


2  tt  Ns  <  r/  > 


^  +  K, 


(9) 


where  r*  =  radius  of  bubbles  measured  directly  in  a  plane.  Equation  (9)  can  be  used  to  calculate  Cs 
which  depends  only  slightly  on  time  even  though  <  r  >  is  a  function  of  time  due  to  bubble  coales¬ 
cence.  No  experimental  data  are  given  in  the  paper  by  deVries,  but  some  references  are  given  by 
Lemlich  (1972). 


Cheng  and  Lemlich  (1983)  studied  errors  which  may  occur  when  the  bubble  size  distribution  in 
a  foam  is  measured  visually  at  a  plane  boundary  of  the  foam.  They  discuss  some  effects  of  statistical 
bias,  bubble  distortion  by  a  solid  contact  plane,  bubble  segregation  at  a  solid  contact  plane,  and 
differences  in  stability  of  bubbles  of  different  sizes.  Planimetric  area  measurements  were  compared 
with  longest  chord  and  longest  median  measurements.  The  planimetric  measurements  were  found -to 
be  the  most  reliable  and  were  in  good  agreement  with  the  true  bubble  size  distribution  for  homoge¬ 
neous  new  foams.  For  heterogeneous  foams,  all  methods  underestimate  the  true  mean  bubble  diame¬ 
ter.  The  main  reason  seems  to  be  that  the  small  bubbles  wedge  themselves  between  the  large  bubbles 
and  the  solid  boundary  plane,  thus  pushing  them  away  from  the  wall.  These  authors  write  Eq.  (7)  as 

7 (r)  ] 


F(r)  = 


The  mean  bubble  radius  for  a  foam  can  be  determined  by  calculating  the  moments  of  the  bubble  size 
distribution  as 


rJ.k  = 


C  rJ  F(r)drJ~k 
J0  rk  F(r)dr 


(11) 
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For  the  biased  distribution  measured  at  the  boundary  surface, 


rs,j  —  1  —  1 


C  rj-lf(r)drj~k 

-U — jjj - - - 

!0  rk~lf (r)dr 


(12) 


The  integrals  are  usually  approximated  by  sums  when  the  mean  radii  are  calculated  from  experimental 
measurements.  The  unbiased  and  biased  mean  radii  are  related  by 


rj,k  =  rsJ-\,k- 1-  (13) 

These  last  four  equations  are  written  using  the  notation  of  Cheng  and  Lemlich  (1983). 

These  authors  generated  foams  with  known  bubble  size  distributions  and  measured  the  mean 
radius  by  observing  the  foam  through  a  glass  boundary.  They  found  that  the  sizes  of  bubbles  observed 
in  the  pool  of  liquid  below  the  foam  corresponded  closely  to  the  true  sizes  of  the  bubbles  in  the  foam. 
That  is, 


r i,o ,p  —  r i,o,/  Sp  ~  fy.  (14) 

Their  data  show  that  as  the  range  of  bubble  sizes  present  becomes  larger,  the  deviation  between  the 
true  mean  bubble  size  and  the  mean  bubble  size  measured  at  the  boundary  plane  also  becomes  larger. 
They  present  their  results  graphically,  and  remark  that  the  graphs  can  be  used  to  predict  the  error 
inherent  to  bubble  size  measurements  made  through  a  boundary  wall.  These  results  are  shown  in  Fig. 

5  where  the  ratio  of  observed  to  true  mean  bubble  radius  is  plotted  as  a  function  of  the  observed  stan¬ 
dard  deviation  of  the  bubble  sizes  made  dimensionless  using  the  observed  mean  bubble  radius.  Figure 

6  shows  the  relationship  between  the  true  standard  deviation  and  the  observed  standard  deviation.  All 
these  results  are  for  newly  formed  foams.  In  these  figures,  the  notation  is  as  follows: 

rs,o,-ij  =  Mean  bubble  radius  in  foam  measured  at  the 
boundary  surface; 

ri,o,i»  =  Mean  bubble  radius  in  the  pool  of  liquid 

below  the  foam; 

f  i,o ,/  =  True  mean  bubble  radius  in  the  foam; 

Sp  =  Standard  deviation  for  bubble  sizes  in  the  pool 

of  liquid  below  the  foam; 

Ss  f  =  Standard  derivation  for  bubble  sizes  measured 

in  foam  measured  at  the  boundary  surface. 

Desai  and  Kumar  (1983)  studied  liquid  hold-up  in  foams  which  were  formed  in  a  polymethyl 
methacrylate  cylinder  by  blowing  air  continuously  through  a  pool  of  liquid  so  that  a  steady-state  foam 
was  established.  The  liquid  hold-up  e  evidently  defined  as  the  volume  fraction  of  liquid  in  the  foam, 
was  found  to  become  smaller  as  distance  from  the  foam-liquid  interface  increased  and  was  found  to 
depend  on  superficial  gas  velocity,  liquid  viscosity,  and  surface  viscosity.  The  value  of  rapidly 
approached  an  asymptotic  value  at  distances  greater  than  approximately  2  cm  from  the  liquid-foam 
interface.  The  asymptotic  value  depended  on  the  gas  velocity  and  the  liquid  properties.  Typical  results 
are  shown  in  Fig.  7.  A  mathematical  model  was  developed  to  predict  liquid  hold-up  profiles  in  the 
steady-state,  semi-batch  foams.  A  unique  feature  of  the  model  is  that  the  liquid  flow  rates  at  the 
foam-liquid  interface  required  for  solving  the  differential  equations  were  calculated  geometrically.  It 
was  shown  that  the  results  of  the  modelling  calculations  are  not  very  sensitive  to  the  boundary  condi¬ 
tions,  so  that  the  geometrical  method  for  estimating  the  boundary  conditions  is  as  effective  as  impos¬ 
ing  experimental  values.  The  results  of  the  modelling  calculations  indicate  that  the  model  predicts  the 
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5  —  Ratio  of  measured  mean  bubble  radius  to  true  mean  bubble  radius 
for  a  newly-formed  foam 
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Fig.  6  —  Relationship  between  true  and  measured  standard  deviations 
for  a  newly-formed  foam 
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Fig.  7  —  Effect  of  superficial  gas  velocity  uc  and  liquid  viscosity 
on  liquid  hold-up  e 
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liquid  hold-up  profiles  for  water  quite  well,  especially  near  the  foam-liquid  interface.  For  higher 
viscosity  liquids,  the  predictions  are  not  as  good. 

Davies  and  Rideal  (1963)  describe  some  simple  methods  for  describing  foams  and  foam  proper¬ 
ties,  such  as 


Expansion  = 


V\ 


(15) 


where  V2  =  volume  of  the  foam; 

Vx  =  original  liquid  volume. 

The  foam  density,  p,  can  be  related  to  the  liquid  density,  pL ,  by 

Vi 

P  =  Pl  TT- 


(16) 


The  foam  viscosity  can  be  up  to  100  iiL  as  measured  in  a  viscometric  flow  device,  because  the  bub¬ 
bles  must  be  deformed,  because  surface  tension  effects  and  surface  viscosity  are  important,  and 
because  the  films  separating  the  bubbles  must  be  detached  from  solid  surfaces  as  the  foam  flows 
through  the  viscometer. 

Princen  (1983)  presents  an  extensive  discussion  of  foam  rheology  and  develops  a  model  relating 
shear  stress  and  shear  rate  for  an  idealized  foam  composed  of  infinitely  long  cylindrical  bubbles.  This 
model  incorporates  elastic  deformation  up  to  a  yield  stress.  For  stresses  greater  than  the  yield  stress, 
the  foam  flows,  and  an  apparent  viscosity  can  be  calculated.  The  model  is  similar  to  the  non- 
Newtonian  Bingham  plastic  model  sometimes  used  to  describe  the  rheology  of  pastes  and  suspensions 
of  solids.  The  yield  stress  and  the  shear  modulus  are  direcdy  proportional  to  the  interfacial  tension 
and  inversely  proportional  to  the  bubble  radius  for  monodispersed  foams.  The  yield  stress  is  a  strong 
function  of  the  volume  fraction  of  gas,  and  the  shear  modulus  is  proportional  to  the  square  root  of  the 
volume  fraction.  The  model  was  first  developed  for  zero  film  thickness  and  zero  contact  angle.  Then 
the  model  was  extended  to  cover  the  more  realistic  cases  of  non-zero  film  thickness  and  contact  angle. 
In  most  cases,  as  the  contact  angle  increases,  yield  stress  increases  and  the  shear  modulus  decreases. 
Yield  stress  and  shear  modulus  are  higher  for  non-zero  film  thickness  than  for  zero  film  thickness. 

An  approximate  measure  of  the  foam- forming  characteristics  of  any  particular  liquid,  the  “foam¬ 
iness,”  can  be  determined  using  a  transparent  tube  fitted  at  the  bottom  with  a  porous  septum.  Gas  is 
blown  upward  through  the  septum  and  through  the  liquid  placed  above  the  septum.  The  steady-state 
volume  of  foam  formed  is  measured  at  constant  volumetric  gas  flow  rates.  The  foaminess  is  defined 
by 


where  V2  =  steady-state  foam  volume; 

u  =  volumetric  gas  flow  rate; 

£  =  foaminess. 

Typical  values  for  £  are  on  the  order  of  5  (clean  seawater)  to  3600  (surfactant  solutions)  seconds 
depending  on  the  liquid  properties. 
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The  stability  of  a  foam  is  a  property  which  is  of  interest  in  many  applications.  A  characteristic 
lifetime  for  a  foam  can  be  defined  by 


‘f-y-i  no  dt 
v  2 


(18) 


where  V2  -  initial  foam  volume; 

V(t)  =  foam  volume  at  any  time  after  formation. 

The  volume  of  the  foam  changes  with  time,  because  liquid  drains  from  the  lamellae,  and  the  bubbles 
break  and  coalesce  due  to  mechanical  shocks  or  vibrations  induced  by  flow  or  surface  tension  effects. 

Haas  and  Johnson  (1968)  have  developed  a  mathematical  model  for  the  drainage  of  liquid  in  a 
stable  foam,  i.e.  one  for  which  the  bubble  size  distribution  is  not  changing.  The  liquid  flows  from  the 
films  separating  the  bubbles  into  the  plateau  border  regions  and  then  through  the  borders  to  the  foam- 
liquid  interface  below  the  foam.  In  the  films,  the  flow  is  driven  by  surface  tension  forces.  In  the  pla¬ 
teau  borders,  the  flow  is  driven  by  gravity.  The  mathematical  model  is  based  on  liquid  phase  mass 
balances  for  the  foam  phase.  For  a  stationary  foam 


e2  _  32^4  7T k2 
LqN  ng  4 


(19) 


where 


Lq 
N(  1 
k 

€pb 


=  volume  fraction  of  liquid  in  the  foam; 

=  superficial  liquid  velocity  for  a  stationary  foam; 

—  e)  =  number  of  plateau  borders  per  unit  cross-sectional  area; 
=  e/e/>B  ; 

=  fractional  cross-sectional  area  for  liquid  flow  in 
plateau  borders. 


The  governing  equation  for  flow  in  the  foam  is 


de  _  d_Lp 
dt  dz 


(20) 


the  position  coordinate  z  is  measured  vertically  upward.  When  Eqs.  (19)  and  (20)  are  combined  to 
eliminate  L0,  the  resulting  differential  equation  can  be  solved  by  the  method  of  separation  of  variables 
to  give  e(t ,z)  =  T(t)Z  (z).  The  solution  is 


e(f  ,z) 


e('o^o)  ~  j  c  To 

32  M 

ng 

v.  j 

irM2 

4 

^  J 

1 

O 

cT0(t  -  t0)  + 

1 

(21) 


where  c  =  integration  constant; 
To  =  T(t  0); 

to  =  initial  value  of  time; 
Zq  =  reference  position. 


Equation  (21)  can  be  combined  with  Eq.  (19)  to  determine  L0.  The  total  amount  of  liquid  drained 
from  the  foam  in  the  time  interval  (t  —  t0)  at  any  position  z  is 
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(22) 


W{t,z)  =  f  L^Lt. 

ro 

At  z  =  Zq,  the  equations  can  be  simplified  to  give 

W(t  ,Zq)  =  L0(t0,zo) 


t  —  tr 


cT0(t  —  t0)  - hi’ 


(23) 


or 


W(t,zQ)  -  -^r -  e2t(t0,z0)  — 

32 M  nNk 2  cT0(t  -  t0)  +  1 


(24) 


The  quantity  cT0  is  evaluated  in  terms  of  the  height  of  the  column  of  foam,  l,  and  the  superficial 
liquid  velocity,  L(r0,z0),  as 


cT  =  4 

0  32 M  xM2 


g2  Lih’Zo)1'2- 


(25) 


The  parameters  l  and  L(r0,z0)  must  be  determined  experimentally.  The  value  of  N  is  related  to  the 
mean  bubble  diameter  by  iV  =  2.5 d~2.  The  experimental  data  reported  by  Haas  and  Johnson  indi¬ 
cate  that  the  model  gives  a  good  description  of  the  drainage  of  liquid  from  a  foam. 


The  foregoing  development  shows  that  the  time  interval  over  which  a  foam  persists  depends  on 
the  viscosity  of  the  liquid.  DeVries  (1972)  reports  lifetimes  of  15  minutes  for  a  foam  formed  from  a 
liquid  having  a  viscosity  /*,  =0.5  poise  but  lifetimes  of  “only  a  few  minutes”  for  foams  formed 
from  water  (jiL  =  0.01  poise).  In  each  case,  VG/VL  =  10.  As  a  foam  ages,  the  bubble  size  distribu¬ 
tion  changes,  and  this  affects  the  foam  stability.  There  are  mechanisms  which  cause  changes  in  bubble 
size  distribution— gas  diffusion  between  bubbles  and  bubble  wall  rupture  and  coalescence.  The  pres¬ 
sure  difference  between  small  and  large  bubbles  is  given  by  Eq.  (2).  This  means  that  there  is  gas 
transport  through  the  bubble  wall  from  the  small  bubble  to  the  large  one.  Eventually,  the  small  bubble 
disappears.  The  rate  of  the  gas  transport  is 


dn 

dt 


PAl2&p  12 


where  P 
A12 

Apl2 


=  permeability  of  the  bubble  wall; 


=  area  of  the  bub 


jle  wall; 


_L_ 

R 


r  =  radius  of  the  small  bubble; 
R  =  radius  of  the  large  bubble; 
n  —  moles  of  gas  transferred. 


The  permeability  can  be  calculated  as 


P 


DS 


9  + 


2D 


(26) 


(27) 
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where  D  =  diffusion  coefficient  of  the  gas  in  the  liquid; 

S  =  solubility  of  the  gas  in  the  liquid; 

9  =  bubble  wall  thickness; 

k  i  =  mass  transfer  coefficient  at  the  gas-liquid  interface 


defined  by 


~  *1  ^i2(Q)  —  C); 


C0  =  equilibrium  concentration  at  the  interface; 
C  =  concentration  in  the  gas  phase. 

There  are  two  limiting  cases. 


(28) 


P 


PS 

9 


.,  2D 

if  — —  «  9 


(29) 


which  applies  if  9  10  4  cm,  and 

p  =  yM  if  ^  »  9.  (30) 

For  soap  films,  9  =  4  x  10  7  cm,  so  that  values  of  k ,  determined  from  flux  measurements  for  soap 

films  using  the  assumption  of  Eq.  (18)  show  that  —  =2  x  10  ~6  cm 

*i 

Diffusion  rate  information  can  be  related  to  bubble  size  distributions.  Let  be  approximately 
the  area  of  the  small  bubble.  The  excess  pressure  in  the  small  bubble  then  is  approximately 

4Pi2  =  (31) 

Combining  Eqs.  (26),  (27),  and  (31)  with  the  ideal  gas  law  shows  that  the  rate  of  shrinkage  of  a  small 
bubble  depends  approximately  on  its  own  radius. 


dr  „ 

-r-=2TP 


(32) 


where  R  =  gas  law  constant; 

T  =  absolute  temperature; 

Pa  =  absolute  pressure. 

Integrating  to  obtain  an  equation  for  r(t), 

r2(0)  -  r\t)  =  rl  -  r2  =  Pt.  (33) 

*  < a 

This  shows  that  the  bubble  surface  area  decreases  as  a  linear  function  of  time.  The  bubble  lifetime  is 
calculated  at  r  =0,  as 
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t(.r0)  = 


-l 


(34) 


4  RTy 

Pa 


r2 

ro- 


DeVnes  (1972)  reports  that  for  one  case,  —  (r2)  *  /(r0).  This  implies  that  d  =  constant  and 

at 

d  ^ 

that  P  —  constant.  For  another  case,  slight  variations  in  —  (r )  with  r0  were  found.  This  could  have 

dt 

been  due  to  variations  in  d  or  P .  If  the  permeability  P ,  is  constant  for  a  particular  foam,  the  change 
in  size  distribution  frequency  function  F(r0),  can  be  calculated  after  selecting  a  method  for  describing 
the  distribution  of  permeated  gas  among  larger  bubbles.  For  example,  the  gas  might  be  distributed 
among  the  larger  bubbles  in  proportion  to  their  surface  areas. 

The  number  density  for  bubbles  is  related  to  the  initial  number  density,  N0,  by 


N, 


r0,t 


\  =  Wo(l  “  ’  F  (r 0)  dr0) 


(35) 


where  F(r0) 

ro,r 

Ndt 


=  frequency  distribution  function  at  t  —  0; 
=  initial  radius  of  a  bubble  with  lifetime  t ; 
—  number  density  at  time  t . 


If  F(r0)  is  known,  then  N it)  can  be  calculated  exactly  knowing  rQ(t)  from  Eq.  (36).  For  a  particular 
foam  formed  mechanically  in  a  high  speed  mixer,  DeVries  (1958)  found  that 


P(r  o)  = 


6a  r0 


(1  +  a  r<5  )4  ’ 

where  a  is  an  adjustable  parameter.  This,  together  with  Eqs.  (34)  and  (35)  gives  the  result 


(36) 


N,  =N0 


1 


1  +  Lt 


,  kd  = 


4  RT 


7  P  a 


(37) 


where  kd  is  in  the  range  10  3  <  kd  <  3.3  x  10  3  seconds. 


-3 


Bubble  wall  rupture  and  subsequent  coalescence  causes  a  spontaneous  decrease  in  surface  area  in 
the  foam  as  time  goes  on.  Wall  rupture  depends  on  wall  thickness  and  on  the  probability  of  rupture. 
Thin  films  are  more  likely  to  rupture  than  thick  ones.  Bubble  wall  lifetime  is  determined  by  the  rate 
of  thinning  and  the  probability  of  rupture.  For  a  bubble  wall  stabilized  by  rigid  monolayers,  flow 
patterns  are  relatively  simple  and  thinning  is  relatively  slow,  where  the 

q2 

(Rate  of  thinning)  —  — . 

HL 


For  a  mobile  bubble  wall,  thinning  is  irregular.  Thin  spots  form  if  the  wall  is  relatively  persistent.  If 
the  wall  is  very  unstable,  it  thins  and  collapses  rapidly  after  reaching  a  critical  thickness  of  about 
10-5  cm.  In  such  a  case 


(Bubble  wall  lifetime)  —  (Rate  of  thinning)  —  \iL . 
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Bubble  wall  rupture  can  be  initiated  by  the  formation  of  a  hole.  Surface  energy  is  converted  to 
kinetic  energy  at  the  growing  hole  rim.  The  velocity  of  hole  rim  propagation  is 


2y 

ur  =  ,  =  constant 


(38) 


For  9  =  10-4  cm,  ur  =  103  cm/s.  This  means  that  rupture  occurs  only  milliseconds  after  initiation 
or  hole  formation.  Thus,  initiation  is  the  rate  determining  step  for  rupture.  Hole  formation  is 
apparendy  necessary  for  rupture  to  occur.  DeVries  (1972)  gives  a  discussion  of  possible  causes  for 
hole  formation.  One  possible  cause  is  that  van  der  Waals  forces  cause  fluctuations  in  film  thickness  to 
propagate.  These  thickness  fluctuations  are  caused  by  thermal  stresses  or  by  some  other  mechanism. 
The  thickness  fluctuations  grow  spontaneously  when  the  wavelength  of  the  fluctuation,  X,  is  greater 
than  some  critical  value.  That  is. 


X  >  ^crir 


2^7 

d^G 

dd2 

L  J 


(39) 


where  G  is  the  total  free  energy  of  interaction  per  unit  surface  area.  Thickness  fluctuations  having 
X  <  \cnt  are  damped. 


Nishioka  et  al.  (1983)  measured  foam  stability  by  observing  the  pressure  in  the  head  space 
above  a  foam.  The  rate  of  increase  of  pressure  was  related  to  the  decrease  in  surface  area.  They 
studied  foams  which  had  different  drainage  characteristics.  They  compared  their  measurements  with 
calculations  made  using  the  mathematical  model  of  Lemlich  (1978).  Typical  results  are  shown  in  Fig. 
8  The  parameter  K  is  a  measure  of  permeability  of  the  liquid  films.  Permeability  increases  as  K 
becomes  larger. 

Djabbarah  and  Wasan  (1985)  studied  the  relationship  between  surface  viscosity  and  foam  stabil¬ 
ity.  They  measured-  the  rupture  time  and  found  that  rupture  time  increased  as  surface  viscosity 
increased.  The  surface  viscosity  was  modified  by  using  surfactant  solutions  containing  differemt  con¬ 
centrations  of  sodium  lauryl  sulfate  and  lauryl  alcohol. 


There  are  other  forces  acting  on  the  bubble  walls  which  can  affect  their  stability.  Forces  associ¬ 
ated  with  an  electrical  double-layer,  viscoelastic  forces,  and  forces  associated  with  surface  tension 
gradients  all  have  an  effect.  There  also  may  be  time  varying  concentrations  of  surfactants  absorbed  on 
the  bubble  wall.  Such  gradients  in  surfactant  concentration  can  cause  surface  tension  gradients  which 
can  either  increase  or  decrease  stability.  If  the  surface  tension  is  lower  at  a  point  on  the  wall  than  it  is 
at  surrounding  points,  surface  energy  can  be  increased  by  extension  of  the  area  of  lower  surface  ten¬ 
sion  and  contraction  of  the  area  of  higher  surface  tension.  This  can  lead  to  rupture  if  the  rate  of  thin¬ 
ning  is  higher  than  the  rate  of  decrease  of  the  surface  tension  gradient.  The  surface  tension  gradient 
would  normally  be  reduced  by  surface  diffusion  of  the  surfactant.  If  this  diffusion  is  too  slow  to 
restore  equilibrium,  then  local  thinning  of  the  bubble  wall  can  lead  to  rupture. 

AIR  ENTRAINMENT 


As  was  mentioned  earlier,  agitation  is  one  method  for  mixing  air  into  a  liquid.  The  turbulent 
flow  in  a  breaking  wave  is  one  type  of  strong  agitation  suitable  for  foam  formation.  This  has  been 
discussed  by  Longuet-Higgins  and  Turner  (1974)  who  developed  a  heuristic  mathematical  description 
of  a  “spilling”  breaker.  They  calculated  approximate  crest-region  shapes  and  discussed  how  the 
length  of  the  whitecap  on  the  breaker  is  related  to  the  energy  balance  in  the  wave.  The  rate  of  air 
entrainment  at  the  base  of  the  whitecap  is  a  topic  left  for  future  work.  More  recently  Griffin  (1984) 
has  discussed  the  overall  aspects  of  wave  breaking,  including  the  formation  of  white  water  and  foam. 
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SURFACE  AREA  (mz) 
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TIME  (hrs) 

Fig.  8  —  The  experimental  decay  of  a  0.4%  sodium  dodecyl  sulfate  foam  compared  with  computed  decays. 
The  value  K  =  10  _4  cm2/hr  represents  the  theoretical  decay  of  a  nondraining  system. 
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Examples  of  the  air  entrainment  and  formation  of  bubbles  during  the  breaking  of  steep  surface  waves 
are  given  by  Su,  Green  and  Bergin  (1984).  A  typical  example  of  a  spilling  breaker  is  given  in  Fig.  9. 
The  region  of  white  water  containing  entrained  air  is  clearly  visible. 

Mixing  air  and  water  by  means  of  plunging  jets  and  by  the  action  of  turbulent  eddies  in  water 
flowing  over  spillways  has  been  discussed  by  several  authors.  Data  for  the  entrained  air  flow  rate  are 
presented  by  Van  De  Sande  and  Smith  (1972,  1973,  1976).  They  directed  a  jet  of  water  at  the  surface 
of  a  pool  of  water  at  such  an  angle  that  the  entrained  air  was  carried  into  a  collection  chamber  con¬ 
nected  to  a  vacuum  system.  The  entrained  air  flow  rate  was  measured  directly,  and  the  data  were 
presented  in  terms  of  the  ratio  of  the  flow  rate,  ctj ,  to  the  jet  flow  rate,  o; .  Their  experimental 
results  can  be  summarized  as  follows.  For  low  velocity  jets  for  which  the  jet  velocity,  Vj,  is  in  the 
range  2m/sec  <  Vj  <  5m/sec,  there  are  two  relationships,  depending  on  the  length  of  the  jet.  These 
are: 


Short  jets:  Lj  <  0.9  LB, 

<t>A  0.027  Vj 

4>j  sin  a 


where  =  320 D„5Vj, 


pl/3 

77172  (SI  units). 
Uj 


Long  jets:  L,  >  0.9  LB, 


<h  =  8.5  x  10~6 
<t>j  <t>j 


+  (1.05  x  10 ~4)  pwVf  (SI  units) 


(40) 


(41) 


where 


h 

Lb 

D: 


a 


D„ 


=  jet  length; 

=  jet  length  at  break-up; 

=  jet  diameter; 

=  angle  between  jet  and  the  liquid  surface; 
=  density  of  the  water; 

=  nozzle  diameter. 


For  high  velocity  water  jets  for  which 


Vj  > 


10  7 

P  A  Dn 


1/2 


(42) 


where  y  =  surface  tension, 
pA  =  density  of  air, 

the  results  are  presented  in  graphical  form  (see  Fig.  2  in  the  1973  paper  by  Van  De  Sande  and 
Smith).  These  results  are  sketched  in  Fig.  10.  The  jet  diameter  is  related  to  the  nozzle  diameter  by 


Dn 

D; 


=  0.085 


X"  —\ 

r  "a 

1/6 

PAVj2Dn 

Pw  VjDn 

=  0.085 

N\Ve  ^  Re 

^  y 

l  ^  J 

1/6 


(43) 


for  NWeNRe  >  7  x  106. 
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Fig.  9  —  A  laboratory  photograph  of  a  deep  water  spilling  breaker; 
from  Ramberg  and  Griffin  (1987) 
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LOG  Vj  (m/s) 

Fig.  10  —  Air  entrainment  by  high  velocity  jets;  from  Van  De  Sande  and  Smith  (1973) 
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The  effect  of  the  angle  of  impact  on  the  air  entrainment  rate  is 


0,4 (sin  o;)I/4  =  constant. 


The  mechanism  for  entrainment  is  described  in  the  1972  paper  by  Van  De  Sande  and  Smith.  Cavities 
are  formed  by  droplet  or  jet  impact,  and  bubbles  are  enclosed.  A  rough-surfaced  jet  entrains  more  air 
than  a  jet  having  a  smooth  surface.  Burgess,  Molloy  and  McCarthy  (1972)  discuss  the  surface  rough¬ 
ness  of  a  jet  and  relate  the  roughness  to  the  amount  of  air  entrained.  Thus, 


,  C  2 

_ D*_ 

<t>j  Dn 


(44) 


where  D*  is  the  diameter  of  the  envelope  which  encloses  all  the  roughness  elements  of  the  jet.  The 
diameter  D*  depends  on  Lj ,  the  length  of  the  jet,  since  the  jet  expands  as  it  moves  away  from  the 
nozzle.  Burgess  et  al.  also  show  how  the  interfacial  area  of  bubbles  entrained  is  related  to  the  jet 
velocity  and  surface  roughness.  These  data  were  determined  by  CC>2  entrainment  and  subsequent 
reaction  with  sodium  hydroxide  solution.  They  calculated  the  interfacial  area  according  to  the  equation 


Rg 


(45) 


where  a 
D 

Ra 

C* 

bo 

k2 


interfacial  area  per  unit  volume  in  the  entrainment  zone; 
diffusion  coefficient; 
reaction  rate; 

equilibrium  CO  2  concentration; 

NaOH  concentration  in  the  liquid; 
reaction  rate  constant; 


and  found  that  in  SI  units 


a  =  (4.5  x  10~3) 


r 

D* 

~K 


Vj  -  (1.75  X  10“2). 


(46) 


The  effects  of  jet  surface  roughness  on  gas  entrainment  are  also  discussed  by  McKeogh  and 
Ervine  (1981).  They  report  the  entrainment  rate  to  be 


<t>A 

+j 


=  1.4 


(47) 


where  e  =  one-half  jet  surface  roughness; 
r  =  mean  jet  radius. 

They  determined  the  size  of  the  region  in  the  liquid  which  was  affected  by  air  entrainment  and  mea¬ 
sured  the  concentration  of  air  in  that  region.  The  concentration  of  air  was  expressed  as 

£,  _  _ 100  x  Volume  of  air _ 

Volume  of  air  +  Volume  of  water 
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The  depth  of  penetration  of  the  entrained  air  was 


Dp  =  2.6  (VjDj)0J. 


The  center  line  air  concentrations  are  given  by 


Q 

r 

max 


7* - ■> 

3 

1  +  3 

X 

Dp  \ 

r  J 

(48) 


(49) 


where  x  =  distance  below  the  surface. 


The  radial  dependence  of  air  concentration  is  given  for  two  regions 


where  rx 


radial  extent  of  the  bubbly  region. 


(50) 


(51) 


Lara  (1979)  has  discussed  the  conditions  necessary  for  air  to  be  entrained  by  a  plunging  jet.  He 
finds  that  whether  or  not  air  entrainment  occurs  depends  on  the  jet  length  to  diameter  ratio  and  on  the 
jet  velocity.  His  results  are  summarized  in  Fig.  11.  The  jet  velocity  is  much  lower  than  the  jet  veloci¬ 
ties  used  in  the  other  work  discussed  here. 


Research  on  air  entrainment  in  water  flowing  over  a  spillway  has  been  reported  by  Gan- 
gadharaiah  et  al  (1970),  by  Keller  and  Rastogi  (1975),  and  by  Wood  (1983).  Their  work  centers  on 
mathematical  modelling  of  the  turbulent  boundary  layer  along  the  spilling  surface.  They  believe  that 
air  entrainment  begins  when  the  outer  edge  of  the  boundary  layer  coincides  with  the  surface  Their  cal¬ 
culations  agree  with  experimental  results  and  observations  made  for  full-size  operating  spillways. 
They  report  the  position  at  which  air  entrainment  begins  as  a  function  of  flow  rate.  Gangadharaiah  et 
al  also  report  the  concentration  of  air  in  aerated  flow  in  a  water  channel,  but  they  do  not  relate  this  to 
the  turbulence  characteristics  of  the  flow. 


Volkart  (1982)  studied  air  entrainment  in  a  partly  filled  pipe  and  related  his  results  to  air 
entrainment  in  open  channel  flow.  The  mechanism  suggested  for  the  air  entrainment  is  that  turbulent 
impulses  perpendicular  to  the  surface  cause  water  droplets  to  be  ejected  into  the  air.  When  these  dro¬ 
plets  fall  back  into  the  water,  they  entrain  air.  In  order  for  water  droplets  to  be  ejected,  the  turbulent 
boundary  layer  must  have  reached  the  surface,  and  the  turbulent  impulse  in  the  transverse  direction 
must  exceed  a  critical  value.  He  postulates  that  this  occurs  somewhat  farther  downstream  than  the 
point  at  which  the  turbulent  boundary  layer  just  reaches  the  surface.  Apparently,  some  extra  distance 
is  required  for  the  turbulence  intensity  to  reach  the  required  value.  The  concentration  of  entrained  air 
in  water  was  measured  using  an  electrical  conductivity  method  based  upon  the  difference  between  the 
conductivity  of  water  and  air.  The  mean  air  concentration  was  correlated  by  the  equation 


C  = 


Qa  +  Qh 


=  1.0  - 


1.0 


(0.02  Ngou  —  6.0)  +1.0 


N, 


BOU 


>  6; 


(52) 
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L/D; 


Fig.  11  —  Conditions  for  air  entrainment  by  jets;  adapted  from  Lara  (1979) 
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where  Qa 
Qw 

Nbou 

Vw 


-  volumetric  flow  of  air  in  the  mixture; 

= volumetric  flow  of  water  in  the  mixture; 

Vw 

=Boussinesq  number  =  — : - 

= water  velocity  with  air  bubble  entrainment 
excluded; 

=mean  hydraulic  radius. 


Equation  (52)  shows  that  air  is  entrained  only  if  NBOU  6.  Volkart  determined  that  there  is  a 
minimum  slope  for  which  NB0U  =  6  and  for  which  air  entrainment  can  be  initiated.  This  slope  is 
given  by 


J A  =  4 


l 


-  log 


2.51? 


+ 


D  'J2gDJa  3.71Z) 


H_ 

D 


=  1.0 


(53) 


where  JA  =  slope; 

ks  =  roughness  height; 
H  =  water  depth. 


The  relationship  giving  friction  losses  for  flow  in  rough  pipes  described  by  Colebrook  (1938)  was 
used  here.  For  depth-to-diameter  ratios  different  from  1.0,  a  correction  factor  is  applied. 


Ja 


~  J  A 


H_ 

D 

V . 


(54) 


The  correction  factor  f(H /D)  is  obtained  graphically.  Typical  values  are 


H_  m 
D  : 


1.0  0.9 


0.8  0.6 


0.5  0.4  0.35 


/-  "N 


V-  J 


1.0  0.94  0.94  0.96 


1.0  1.05 


1.08. 


The  initiation  of  air  entrainment  depends  on  pipe  diameter,  roughness  height,  and  depth  to  diam¬ 
eter  ratio.  The  results  that  air  entrainment  begins  at  lower  values  of  JA  for  large  smooth-walled  pipes 
than  for  small  rough-walled  pipes.  The  author  says  this  is  evidence  that  turbulence  is  a  more  impor¬ 
tant  factor  than  wall  roughness  when  it  comes  to  describing  air  entrainment.  The  critical  velocity  for 
air  entrainment  to  occur  is 


K  .crit  .  (55) 

The  mean  velocity  of  the  entire  air-water  mixture  was  found  to  be  smaller  than  the  mean  velocity  for 
the  water  only,  or 


(56) 
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The  air  concentrations  for  flow  in  partly  filled  pipes  are  smaller  than  those  for  open  channel  flow  at 
the  same  conditions  of  slope,  roughness  height  and  depth  to  diameter  ratio. 

MARINE  HYDRODYNAMICS  APPLICATIONS 

Wave  Breaking  There  are  many  causes  for  the  breaking  of  ocean  waves.  These  include  the  rela¬ 
tive  motion  between  the  water  and,  say,  a  ship  in  a  seaway  or  a  cylindrical  obstacle  in  a  wavefield; 
wind  blowing  over  the  water;  superposition  of  and  nonlinear  interactions  between  wave  components; 
concentration  of  wave  energy  by  refraction;  and  the  shoaling  of  waves.  These  various  causes  result  in 
two  predominant  types  of  breaking  waves:  the  plunging  breaker,  in  which  the  wave  crest  curls  for¬ 
ward  and  plunges  into  the  slope  of  the  wave  at  some  distance  away  from  the  crest;  and  spilling  break¬ 
ers,  in  which  the  broken  region  tends  to  develop  more  gently  from  an  instability  at  the  crest  and  often 
forms  a  quasi-steady  whitecap  on  the  forward  face  of  the  wave.  Foam  and  bubbles  from  the  entrained 
air  often  play  an  important  and  visible  role  in  the  processes  of  spilling  and  plunging. 

A  third  type,  surging,  sometimes  develops  as  waves  are  incident  upon  a  sloping  beach.  In  this 
case,  if  the  slope  is  very  steep  or  the  wave  steepness  is  very  small,  the  waves  do  not  actually  break 
but  surge  up  and  down  to  form  a  standing  wave  system  with  little  or  no  air  entrainment.  A  fourth 
case,  collapsing,  is  considered  to  be  a  special  limiting  case  of  the  plunging  breaker.  Collapsing  occurs 
when  the  crest  remains  unbroken,  but  the  lower  part  of  the  front  face  of  the  wave  steepens  and  falls 
and  then  forms  an  irregular  region  of  turbulent  water.  Galvin  (1968)  has  provided  a  historical 
recently,  the  breaking  of  waves  on  beaches  has  been  discussed  in  some  detail  by  Peregrine  (1983). 

The  breaking  of  waves,  with  particular  emphasis  on  the  processes  which  occur  in  deep  water, 
also  has  been  discussed  by  Griffin  (1984).  From  numerous  small-scale  experiments  it  is  generally 
expected  that  air  entrainment  is  not  a  determining  factor  in  understanding  or  predicting  the  dynamics 
of  breakers  (D.H.  Peregrine,  1984;  private  communication).  However,  air  entrainment  and  the  ensu¬ 
ing  formation  of  bubbles  and  foam  do  give  a  measure  of  the  intensity  and  extent  of  the  turbulence 
which  is  produced  by  spilling  and  plunging  breakers.  Two  examples  are  shown  in  Figs.  12  and  13.  A 
laboratory  simulation  of  the  air  entrainment  and  bubble  formation  processes  during  wave  breaking  is 
shown  in  Fig.  12.  The  formation  of  a  bubble  cloud  beneath  the  breaker  is  clearly  visible,  and  it  is 
evident  that  even  this  nominally  simple  laboratory-scale  flow  is  very  complex.  The  breaking  of  a 
steady  surface  wave  and  the  ensuing  formation  downstream  of  a  turbulent  white  water  wake  region  is 
shown  in  Fig.  13.  The  wave  train  was  generated  by  towing  a  foil-shaped  body  beneath  the  free  sur¬ 
face.  Again  the  flow  field  is  very  complex  for  even  this  simple  configuration.  Further  progress  toward 
modeling  and  understanding  the  turbulent  regions  of  foam  and  bubbles  in  breaking  waves  is  likely  to 
be  empirical  for  some  time  to  come,  except  for  simple  models  of  bubble  clouds  as  described  in  the 
next  section. 

The  extent  of  surface  coverage  by  foam  and  bubbles  is  important  in  characterizing  and  interpret¬ 
ing  the  return  from  remote  sensing  devices  such  as  microwave  radiometers  and  microwave  radars 
which  are  deployed  above  the  sea  surface.  The  radiation  from  the  ocean  surface  that  is  measured  by 
microwave  radiometry  is  a  sum  of  the  thermal  emission  of  the  water  and  the  reflection  of  the  sky’s 
radiation  from  the  sea  surface.  Foam  greatly  increases  the  microwave  emissivity  of  the  water  surface 
(Droppleman,  1970;  Stogryn,  1972;  Koepke,  1986)  by  altering  the  reflectivity  (emissivity  =  1  - 
reflectivity)  due  to  absorption  and  reflection  from  the  upper  air-foam  and  the  lower  foam-water  inter¬ 
faces.  The  foam-induced  increase  in  emissivity  is  primarily  responsible  for  the  large  increase  in  the 
microwave  brightness  temperature  of  the  sea. 

The  presence  of  surface  foam  can  contribute  significant  backscatter  to  microwave  sensors 
(Alpers,  Ross  and  Rufenach,  1981;  Keller,  Plant  and  Valenzuela,  1986).  In  addition  to  the  return 
from  pure  Bragg  scattering,  there  are  contributions  from  specular  scattering  by  the  wave/turbulent 
fluid  interfaces  normal  to  the  incidence  angle  of  the  radar  and  volume  scattering  when  the  foam  is 
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Fig. 


12  —  Sequence  of  air  entrainment  and  bubble  formation  in  a  breaking  wave. 
The  wave  moves  from  left  to  right;  from  Su,  Green  and  Bergin  (1984). 
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Fig.  13  —  A  steady  breaking  surface  wave  and  the  ensuing  turbulent  wake  downsteam 
from  a  submerged,  towed  hydrofoil.  Photo  courtesy  of  J.H.  Duncan. 


well  aerated.  Significant  penetration  of  the  radiation  into  the  foam  takes  place  when  the  dielectric  pro¬ 
perties  of  the  foam  approach  those  of  air.  More  aspects  related  to  whitecaps  and  white  water  genera¬ 
tion  by  breaking  surface  waves  are  discussed  in  the  proceedings  of  a  recent  whitecap  workshop 
(Monahan  and  MacNiocall,  1986). 

Surface  Ship  Wakes.  Figure  14  illustrates  the  major  sources  of  surface  foam  and  bubble  produc¬ 
tion  from  displacement-type  vessels  and  defines  some  of  the  relevant  terminology.  The  surface  wake 
can  be  divided  into  three  regions.  These  are  the  near  wake,  the  far  wake,  and  the  Kelvin  wave  sys¬ 
tem. 


The  Kelvin  wave  system  19.5  boundary  (cusp  lines)  and  35.3  cusp-crest  tangent  lines  are  illus¬ 
trated  in  the  figure.  The  apex  of  the  boundary  lines  is  located  about  one  ship  length  ahead  of  the  bow 
(Newman,  1970).  The  transverse  and  divergent  wave  crests  are  visible  for  many  ship  lengths  astern 
and  to  either  side  of  the  ship’s  path.  Surface  foam  and  bubbles  can  be  produced  when  the  divergent 
waves  in  the  Kelvin  wave  system  break  at  wave  steepnesses  ak ,  where  a  is  the  wave  amplitude  and  k 
-  2 ir/L  is  the  wavenumber,  greater  than  ak  —  0.34  (Ramberg  and  Griffin,  1987).  Energy  input  into 
the  Kelvin  wave  system  by  a  vessel’s  movement  throught  the  water  and/or  the  interaction  of  the  Kel¬ 
vin  wave  and  ambient  wave  fields  also  can  cause  the  wave  steepness  to  exceed  ak  =  0.34.  Figure  1 
shows  breaking  (spilling  type)  divergent  waves  almost  directly  perpendicular  to  the  midsection  on  the 
port  and  starboard  sides  of  the  ship. 

There  are  three  major  sources  of  foam  and  bubble  production  in  the  near  wake  region.  The  bow 
wave  that  is  generated  by  the  ship’s  motion  will  break  if  its  steepness  is  greater  than  about  ak  = 
0.34.  The  turbulent  region  adjacent  to  the  ship’s  hull  produces  foam  and  bubbles  because  of  the  fric¬ 
tional  drag  forces  on  the  surface  of  the  hull.  The  initial  spreading  region  (ISR)  of  the  propeller  wake 
is  a  region  of  high  angular  divergence  of  foamy  water  directly  aft  of  the  ship’s  stem.  It  is  generally 
outlined  by  what  appears  to  be  a  spilling-type  breaking  wave. 

At  greater  distances  downstream  in  the  far  wake,  the  surface  foam  and  bubbles  decay.  The  per¬ 
sistence  of  the  foam  depends  on  the  time  it  takes  for  the  bubbles  to  break  after  they  reach  the  surface. 
The  major  factors  that  increase  the  stability  of  a  bubble  at  the  surface  are  increasing  salinity  (Miyake 
and  Abe,  1948;  Peltzer  and  Griffin,  1987;  Monahan  and  Zeitlow,  1969;  Scott,  1975),  decreasing 
water  temperature  (Miyake  and  Abe,  1948)  increased  surface  viscosity  (Kitchener  and  Cooper,  1959), 
and  the  presence  of  organic  surface-active  materials  which  reduce  the  surface  tension  and  modify  the 
surface  rheology  (Adamson,  1976;  Akers,  1976;  Bikerman,  1973).  A  recent  photographic  analysis 
(Peltzer,  1984)  revealed  that  the  length  of  this  foamy  white-water  region  (Lww)  can  be  described 
equally  well  by  the  equations 

Avw  -  65.5  V^n-0-5 


and  =  53.5  /(gL^)05]1'40 

where  Vs  (m/s)  is  the  ship  speed,  LwL  is  the  waterline  length  (m),  n  is  the  propeller  revolutions  per 
second  and  g  is  the  gravitational  constant  (m/sz). 

A  typical  wake  profile  taken  of  a  surface  vessel  at  ultrasonic  frequency  is  shown  in  Fig.  15, 
which  has  been  adapted  from  a  Naval  Ocean  Systems  Center  (NOSC)  study.  The  wake  begins  at  the 
propellor  and  then  increases  to  a  depth  of  two  to  three  times  the  ship’s  draft.  According  to  the  study, 
the  rapid  wake  depth  increase  is  accompanied  by  very  high  bubble  density  during  the  first  10  to  20 
seconds.  This  probably  represents  the  initial  spreading  region  of  the  wake.  It  should  be  noted  that  the 
wake  width  increase  at  these  early  times  is  much  more  rapid  than  the  estimate  given  in  the  figure. 
The  bubble  density  decreases  at  times  between  20  and  200  seconds,  the  wake  depth  remains  relatively 
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Fig.  15  —  A  general  surface  wake  profile  measured  acoustically  at  ultrasonic  frequency 


constant  at  two  to  three  times  the  ship’s  draft,  and  the  wake  width  is  somewhat  greater  than  the  ship’s 
beam.  At  times  beyond  200  seconds  the  depth  of  the  wake  gradually  decreases  and  the  wake  width 
spreads  to  about  twice  the  ship’s  beam.  The  time  period  between  20  and  500  seconds  probably 
represents  the  far  wake.  For  a  typical  125  m  long  ship  moving  at  5.2  m/s,  a  time  of  500  seconds 
represents  a  downstream  distance  of  11  ship  lengths. 

Modelling  of  the  downstream  surface  ship  wake  will  depend  on  the  availability  of  appropriate 
initial  conditions.  One  logical  place  to  specify  these  conditions  is  in  a  vertical  plane  at  the  end  of  the 
initial  spreading  region.  Then  the  spatial  and  temporal  evolution  of  the  far  wake  region,  in  principle, 
can  be  computed.  However,  the  computations  are  strongly  dependent  on  good  initial  conditions  which 
in  turn  depend  upon  the  accurate  characterization  of  the  several  sources  of  foam  and  bubble  genera¬ 
tion  upstream  of  the  far  wake. 

The  residence  time  of  a  single  or  multilayer  foam  at  the  surface  and  the  bubbles  in  the  underly¬ 
ing  water  column  depends  on  the  stability  of  the  bubbles  in  those  layers.  Some  of  the  major  factors 
that  influence  the  stability  of  bubbles  include  the  water  temperature,  the  air-sea  temperature  differ¬ 
ence,  the  surface  viscosity,  the  presence  of  organic  surface-active  materials,  water  cleanliness,  and  the 
presence  of  biological  organisms  and  particulate  matter.  Other  important  factors  include  the  salt  con¬ 
tent  or  salinity  of  the  water,  the  size  of  the  bubbles,  the  humidity  and  velocity  of  the  air  above  the 
foam  layer  and  the  thickness  of  the  liquid  lamallae.  All  of  these  factors  must  be  considered  when  the 
initial  conditions  are  specified. 

It  will  be  necessary  further  to  specify  the  relative  proportions  of  foam  and  the  “normal”  two- 
phase  mixture  of  bubbles  and  sea  water.  The  foamy  region  is  confined  to  a  relatively  thin  vertical 
layer  just  at  and  below  the  surface.  This  is  a  highly  agitated  and  aerated  region  with  a  high  relative 
content  of  air  to  water.  It  may  be  possible  to  assume  that  the  foam  density  is  constant  in  the  region. 
Then  below  this  is  the  bubble  wake  which  extends  vertically  downward  over  a  greater  distance.  The 
sparse  knowledge  which  is  available  concerning  the  vertical  distance  over  which  the  wake  extends  is 
summarized  in  Fig.  15.  Typically,  acoustic  measurements  of  the  vertical  wake  thickness  have  dealt 
primarily  with  the  bubbly  region.  As  noted  previously  the  thickness  of  this  region  is  approximately 
two  ship  drafts  in  depth.  The  vertical  extent  of  the  overall  turbulent  ship  wake  probably  is  greater 
still.  Models  such  as  the  ones  described  by  Thorpe  (1982)  and  Miner  et  al  (1986)  for  the  motions  of 
individual  bubbles  and  clouds  of  bubbles  in  seawater  may  provide  a  first  step  toward  a  more  complete 
understanding  of  the  bubble-dominated  portion  of  the  turbulent  wake. 
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